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[i] Airborne Multi-AXis Differential Optical Absoiption Spectroscopy (AMAX-DOAS) 
measurements of N0 2 tropospheric vertical columns were performed over California for 
two months in summer 2010. The observations are compared to the NASA Ozone 
Monitoring Instrument (OMI) tropospheric vertical columns (data product v2.1) in two 
ways: (1) Median data were compared for the whole time period for selected boxes, and the 
agreement was found to be fair (R = 0.97, slope = 1.4 ±0.1, N= 10). (2) A comparison was 
performed on the mean of coincident AMAX-DOAS measurements within the area of the 
corresponding OMI pixels with the tropospheric NASA OMI N0 2 assigned to that pixel. 

The effects of different data filters were assessed. Excellent agreement and a strong 
correlation (R = 0.85, slope = 1.05 ±0.09, N = 56) was found for (2) when the data were 
filtered to eliminate large pixels near the edge of the OMI orbit, the cloud radiance fraction 
was <50%, the OMI oveipass occurred within 2 h of the AMAX-DOAS measurements, the 
flight altitude was >2 km, and a representative sample of the footprint was taken by the 
AMAX-DOAS instrument. The AMAX-DOAS and OMI data sets both show a reduction of 
N0 2 tropospheric columns on weekends by 38 ±24% and 33 ± 1 1%, respectively. The 
assumptions in the tropospheric satellite air mass factor simulations were tested using 
independent measurements of surface albedo, aerosol extinction, and N0 2 profiles for Los 
Angeles for July 2010 indicating an uncertainty of 12%. 
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DOAS measurements over California: Testing the NASA OMI tropospheric N0 2 product, J. Geophvs. Res. Atmos., 118, 
7400-7413, doi: 10. 1002/jgrd.50550. 


1. Introduction 

[ 2 ] Nitrogen dioxide (N0 2 ) is a regulated air pollutant, and 
monitoring is needed to assess the effects of reduction 
measures. Satellites provide the opportunity for global 
measurements of tropospheric vertical column densities 
(tropVCDs) on the basis of nadir measurements of ultraviolet 
(UV) and visible scattered sunlight. There is/was a number of 
satellite instruments available for measuring N0 2 tropVCDs 
starting in the mid-90s of the last century: GOME on ERS-2 
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[Burrows et al. , 1999], SCIAMACHY on ENVISAT 
[Bovensmann et a!., 1999], Ozone Monitoring Instrument 
(OMI) on AURA [Levelt et al., 2006], and GOME-2 on 
MetOp-A [Munro et al., 2006], Although satellites today 
only have a relatively coarse resolution of ca. 1 0 to 1 00 s of 
kilometers and usually only provide a global picture every 
1 - 6 days, those measurements are useful tools for trend 
analyses [e.g., Richter et al., 2005; Lamsal et al., 2011], to 
infer NO x emissions by inverse modeling [e.g.„ Miyazaki 
et al., 2012; Lin, 2012], as assimilation data for air quality 
forecasting [e.g., Wang et al., 2011; Petritoli et al., 2011], 
and to deduce N0 2 surface concentrations [e.g., Lamsal 
et al., 2008; Lee et al., 2011], However, tropVCDs from 
different research groups do not always agree [e.g., Lamsal 
et al., 2010; Bucsela et al., 2008; Boersma et al., 2008]. 
Differences result mainly from differences in the assump- 
tions in the radiative transfer calculation and in the way the 
stratospheric N0 2 is removed from the total column density. 
The uncertainty of satellite retrievals for polluted environ- 
ments is dominated by the uncertainty in the radiative 
transfer calculations and originates from the necessary 
assumptions made to describe the surface albedo, the N0 2 
and aerosol extinction a priori profile shapes as well as insuf- 
ficient information about clouds [Boersma et al., 2004; 
Richter and Burrows, 2002]. Therefore, regular validation 
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Figure 1 . BAMFs for satellite and aircraft geometries for a 
SZA of 17° (time of OMI/ AURA overpass over North 
America). Aircraft BAMF was calculated for aircraft 
altitudes of 3 km and 5 km. The BAMFs have been calculated 
for 435 nm for OMI and otherwise for 446 nm (center wave- 
lengths of the NO 2 fitting windows). Aerosol optical depths 
are 0.2 and 0.5 and have a negligible effect on the nadir 
observations. The nonzero zenith BAMF below the aircraft 
highlights the need for choosing a reference spectrum with 
low trace gas concentration. 

under a variety of conditions and also spanning the entire 
area of a satellite footprint is important. 

[3] We present an NO 2 validation study using the 
University of Colorado Airborne Multi-Axis Differential 
Optical Absorption Spectroscopy (CU-AMAX-DOAS) 
instrument, which is optimized to obtain tropVCDs of trace 
gases as well as detailed profile information for trace gases 
and aerosol extinction when flying so-called low approaches 
[Baidar et al., 2013]. Werecorded data in California from 19 
May - 19 July 2010 flying on a NOAA Twin Otter aircraft 
(http://www.aoc.noaa.gov/aircraft_otter.htm) and compared 
the data to the NASA OMI tropospheric N0 2 product. 
The sampling area contained pollution hotspots like the 
Los Angeles Metropolitan Area as well as background 
conditions, e.g., over the High Desert. These measurements 
were part of the California Research at the Nexus of Air 
Quality and Climate Change (CalNex, campaign overview 
by Ryerson et al. [2013]) and Carbonaceous Aerosols and 
Radiative Effects Study (CARES, campaign overview by 
Zaveri et al. [2012]) campaigns. We also deployed a set of 
upward and downward pointing radiometers to measure 
surface albedo. Hence, the key parameters for the solar 
radiative transfer calculation, i.e., surface albedo, aerosol 
extinction profiles, and trace gas profiles, were measured 
quasi-simultaneously with the NCB tropVCD. Since the 


AMAX-DOAS instrument measures in a similar wavelength 
range as the UV/visible satellite spectrometers, the same 
absorbers can be retrieved which makes this technique 
applicable to the validation of other species as well, e.g., 
glyoxal or formaldehyde. 

[4] Satellite instruments and the CU-AMAX-DOAS 
instrument measure with a similar viewing geometry, i.e., 
nadir or near nadir. Hence, they exhibit similar vertical sensi- 
tivities. Figure 1 compares box air mass factors (BAMFs; 
BAMFs are equivalent to weighting functions for optically 
thin absorbers, and they express the average light path 
enhancement with respect to the vertical thickness of a layer) 
for a typical scenario during CalNex/CARES at the solar 
zenith angle (SZA) of the OMI overpass. AMAX-DOAS 
measurements are on average about 50 - 36% more sensitive 
than OMI to the lowest 1 km in the atmosphere for an aircraft 
altitude of 3 - 5 km, respectively, but the shapes of the 
BAMFs are very similar. The aerosol extinction in these 
examples is rather at the extremes: AERONET measure- 
ments at CalTech, Pasadena indicate a mean aerosol optical 
depth (AOD) of 0.19 ±0.08 and a median AOD of 0.1 8 with 
0.14 being the 25th percentile and 0.24 being the 75th 
percentile at 440 nm over the time period of 1 9 May - 1 9 
July 2010. As can be seen from Figure 1, the differences in 
the BAMFS for the two different AODs are negligible for 
the nadir observations (compare top and bottom panel). 

[5] Tropospheric NO 2 satellite validation has been 
performed with airborne or ground-based chemilumines- 
cence NO x analyzers, but also with airborne or ground-based 
remote-sensing instruments. 

[6] When comparing ground-based concentrations with 
column densities measured from space, assumptions of the 
profile shape have to be made either by using models 
[ Petritoli et al., 2004; Ordonez et al., 2006; Blond et al., 
2007; Lamsal et al., 2010] or mixing depth from climatol- 
ogies [Boersma et al., 2009] or as done by Schaub et al. 
[2006], using NO 2 measurements from different altitudes in 
the Alps to construct a profile. The chemiluminescence 
NO x analyzers applied in the above studies use a heated 
molybdenum oxide surface to reduce N0 2 to detectable 
NO. These analyzers are known to suffer from interferences: 
species such as nitric acid or peroxyacetyl nitrate also convert 
to NO [e.g., Steinbacher et al., 2007; Dunlea et al., 2007], 
and correction factors have to be applied [e.g., Ordonez 
et al., 2006], Discrepancies between the different measure- 
ments can partly be attributed to the nonrepresentativeness 
of point measurements [e.g., Ordonez et al., 2006; Schaub 
et al., 2006; Celarier et al., 2008; Boersma et al., 2009]. 
Hence, ground-based in situ monitor validation is best 
performed for background sites [e.g., Lamsal et al., 2010]. 
Especially at urban hotspots, the traditional and localized 
ground-based validation attempt fails [e.g., Lamsal 
et al., 2010]. 

[7] In situ measurements of NCB on an aircraft make 
profile measurements possible when the altitude of the 
aircraft is changed. Usually, spiral flight patterns are 
performed for this. The concentrations can then be integrated 
to yield partial column densities. However, missing layers 
have to be interpolated and the layers above maximum flight 
altitude or below the minimum extrapolated to obtain a 
tropVCD to compare with satellite products [Heland et al., 
2002; Ladstatter-Weifienmaver et al., 2003; Martin et al., 
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2004; Bucsela et al., 2008; Boersma et al., 2008; Hains et al., 
2010], This approach is prone to introduce errors for the 
bottom layer where most of the NO 2 resides and an aircraft 
is not allowed to fly, at least not if not taking off or landing 
at an airport. 

[s] Ground-based multi-axis (MAX) DOAS measure- 
ments have been used for comparison with satellites since 
those instruments not only provide tropVCDs, but also 
surface concentrations for comparison with NO x analyzers 
[e.g., Brinksma et al., 2008], Measurements usually agree 
within the limits of the uncertainties. But the quality of agree- 
ment is variable depending on satellite product, but also 
season and location. MAX-DOAS tropVCDs at highly 
polluted areas seem to indicate that also those measurements, 
similar to the ground-based in situ sensors, are more sensitive 
to local emissions and yield higher columns than the satellite 
(Figure 1) [e.g., Irie et al., 2008; Brinksma et al., 2008; Ma 
et al., 2013; Me et al., 2012], Other remote-sensing 
ground-based techniques that have been used for satellite 
validation include Lidar [Volten et al., 2009], zenith-sky 
DOAS measurements [Petritoli et al., 2004; Chen et al., 
2009], and direct-sun DOAS [Wenig et al., 2008]. 

[ 9 ] Although airborne DOAS NO 2 measurements have 
been performed for several years [e.g., Melamed et al., 
2003; Wang et al., 2005; Merlaud et al., 2011], only in three 
studies, N0 2 tropVCDs were compared to satellite measure- 
ments but with very limited data from only one day of 
measurements each [Heue et al., 2005, 2008; Bruns et al., 
2006; Walter et al., 2012], Knowledge of the horizontal 
variability of tropVCDs can also be obtained by ground- 
based mobile MAX-DOAS, and those tropVCDs were 
compared to OMI but again only for a few days of measure- 
ments [Wagner et al., 2010; Shaiganfar et al., 201 1]. 

[ 10 ] The advantages of using the CU-AMAX-DOAS mea- 
surements recorded during summer 2010 over California for 
NASA OMI N0 2 tropVCD validation are the vast amount of 
data collected over a time period of two months (overall 
206 h) spanning an area covering Los Angeles and surround- 
ings (and to a lesser extent Sacramento and surroundings). 

2. Instrumentation and Measurements 

[ 11 ] The CU-AMAX-DOAS measurements used for this 
validation study were recorded during the CalNex and 
CARES campaigns, both took place in California in summer 
2010. The CU-AMAX-DOAS instrument and albedo sensors 
were mounted to the NOAA Twin Otter remote-sensing 
aircraft which was based out of Ontario International 
Airport during CalNex from 19 May - 15 June and 30 June 
- 19 July and for a shorter time period (16 - 29 June), 
McClellan Airport in Sacramento during CARES. The 
majority of the measurements were perfonned over the 
Greater Los Angeles Area and its relatively unpolluted 
surroundings. Hence, measurements of the study presented 
here provide a large range of NO x conditions from back- 
ground to highly polluted. 

2.1. CU-AMAX-DOAS Instrument and N0 2 
tropVCD Retrievals 

[ 12 ] The instrument and its performance during CalNex 
and CARES are described in detail in Baidar et al. [2013]. 
Briefly, a telescope pylon is window-plate-mounted pointing 


parallel to the flight direction. A motor-mounted prism in the 
telescope permits scanning in forward, upward, and down- 
ward directions under discrete elevation angles with a field 
of view (FOV) of 0.3° x 5.9°. Independently of the move- 
ment of the aircraft, the angle is held constant with respect 
to horizontal by means of an angle sensor coupled to a motor 
adjusted with controller times of~20Hz and an overall 1 
sigma angle accuracy of 0.35°. The collected light is 
transmitted by a split glass fiber bundle into two tempera- 
ture-stabilized spectrometer/CCD detector systems. In this 
study, spectra only from one of the systems are used with a 
wavelength range of 330 - 470 nm and a resolution of 
0.7 nm full-width half-maximum (FWHM). The spectra are 
analyzed with the well-known DOAS technique [e.g., Platt 
and Stutz, 2008]. Here, the WinDOAS software [Fayt and 
Van Roozendael, 2001] was used to retrieve NO 2 slant 
column densities (SCDs) in a wavelength range from 433 
to 460 nm applying two NO 2 absorption cross sections at 
220 K and 294 K [Vandaele et al., 2002], and further includ- 
ing cross sections for the oxygen dimer [Greenblatt et al., 
1990], ozone [Bogumil et al., 2001], glyoxal [Volkamer 
et al. 2005], and water vapor [Rothman et al., 2005], A ring 
cross section was calculated with the MFC software 
[Corner et al., 1993] from the zenith reference spectrum 
and taken into account as a pseudo-absorber in the DOAS 
fit. A polynomial of the order of three was used and a 
straylight correction (see Fayt and Van Roozendael, 2001) 
applied. All spectra of one flight were analyzed in compar- 
ison to one common zenith reference spectrum recorded 
during that particular flight in a clean environment and fly- 
ing at relatively high altitudes, i.e., ca. 3 - 5 km. The results 
are trace gas differential SCDs (dSCDs), i.e., the integrated 
absorber density along the average photon path from the sun 
to the telescope, and differential with regards to the amount 
of NO 2 present in the reference spectrum. 

[ 13 ] In the following, only nadir and zenith spectra are 
investigated, for other geometries see Baidar et al. [2013]. 
Stratospheric NO 2 concentrations are relatively stable but 
undergo photochemical reactions resulting in a distinct 
diurnal pattern with higher N0 2 towards higher SZAs. In 
addition, the dSCDs also increase towards higher SZAs 
due to the longer light path through the stratosphere. 
Therefore, for some early or late measurements, fast 
changes in the stratospheric N0 2 over the course of a flight 
occurred. This change is a smooth function with time of day 
for a certain location. However, changing the location of the 
aircraft also changes the local SZA. A polynomial of the 
order of three was fitted through all the zenith dSCDs 
(the frequency for zenith measurements varied, but was up 
to every 30 s towards the latter half of the campaign, see 
also below) of such a flight with aircraft altitudes of 
1800 m above sea level (ASL) or more. This polynomial 
was then subtracted from the nadir dSCDs which yields 
an effective tropospheric SCD (tropSCD). The polynomial 
approach was chosen instead of a linear interpolation of 
the zenith dSCDs in order to better account for slight 
variations in the possible residual tropospheric NO 2 origi- 
nating from possible tropospheric NCL above the flight al- 
titude as well as changes in stratospheric NO 2 due to the 
slight changes in the local SZA. This leads to a relatively 
high uncertainty of 0.5 - 2.0 x 10 15 molecule cm~ 2 in the 
stratospheric correction when the polynomial had to be 
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applied to the measurements. For many flights, the nadir 
dSCDs can be directly interpreted as a tropSCD assuming 
that the zenith reference spectrum, obtained at altitudes 
3 km, only contains stratospheric N0 2 and that the change 
in stratospheric NO 2 is negligible during a 4h flight. Of 41 
flights, 18 indicated a changing stratospheric NO 2 load 
during parts of the flight, and those were corrected with 
a polynomial. 

[ 14 ] The original or corrected tropSCD is then converted to 
a tropVCD with a geometric air mass factor ( geoAMF ): 


tropVCD 


tropSCD 

geoAMF 


tropSCD 
sscSZA + 1 


( 1 ) 


[ 15 ] For this geoAMF, the assumption is made that the 
photons reach the telescope after being reflected from the 
Earth's surface in the nadir point. Hence, the slant column 
in comparison to the vertical column is weighted by the 
secant function of the SZA plus unity for the part of the 
lightpath from the nadir point into the telescope. The aircraft 
height and hence the integration height is implicitly included 
in this equation by applying the trigonometric function of a 
right-angled triangle. After extensive radiative transfer calcu- 
lations simulating a range of possible scenarios for NO 2 
mixing ratios and mixing heights, optical aerosol parameters, 
AOD, and surface albedo, the geometric approximation was 
found to be an adequate approximation and chosen here to 
be independent of climatology or model NO 2 profile data 
which would be needed as a priori in radiative transfer simu- 
lated AMFs. Baidar et al. [2013] have shown that this 
geometric conversion yields reliable results within 10% 
comparing with ground-based MAX-DOAS instruments. 
The error due to the use of geoAMF depends on the SZA 
and is accurate to within 6% for SZA <20°, 10% for 
SZA <50°, and 20% for SZA <60° when compared to 
radiative transfer simulation performed with the McArtim 
code [Deutschmann et al., 2011] for the predominant condi- 
tions encountered during the CalNex/CARES campaigns 
[ Baidar et al., 2013], Most flights were conducted at 
SZA <60°. The uncertainty on the N0 2 cross section is ca. 
5% depending on the temperature [ Vandaele et al., 2002], 
For the nadir dSCD, the 1 sigma DOAS fitting error is ca. 
3% and the limit of detection ca. 3.2 x 10 15 molecule cm" 2 
for 2 s integration time. In addition, the uncertainty in the 
geoAMF (see above) and the uncertainty of 0.5 - 2.0 x 10 15 
molecule cm" 2 in the stratospheric correction have to be 
taken into account. An overall estimate for the uncertainty 
of the tropVCD is 10 - 20% depending on the SZA. 
Please note that the overpass of the AURA satellite occurs 
at about 17° SZA over California and therefore at a time 
closest to the minimum uncertainty of the AMAX- 
DOAS measurements. 

[ 16 ] The CU-AMAX-DOAS data acquisition software was 
improved over the course of the two months resulting in 
integration times of 2 s for the latter half of the campaigns. 
The scanning sequence when flying at altitudes above 2 km 
above ground level (AGL) included the angles 90°, —2°, 
—5°, -10°, -20°, -90°, 0°, 2°, 5°, 10°, and 20° (positive 
angles upward and negative downward), with a larger 
number of observations on the downward looking angles. 
Typically, nadir spectra were recorded every 12 - 15 s. 


With a typical aircraft speed of 65 ms -1 , this corresponds 
to a horizontal translation of about 900 m. The FOV of the 
telescope gives a nadir footprint of~20 m across and~550 m 
along the track while flying at 4 km altitude for 2 s 
integration times. 

2.2. NASA OMI Instrument and N0 2 
tropVCD Retrievals 

[ 17 ] OMI flies on the NASA Aura satellite which was 
launched into a near-polar sun-synchronous orbit in 2004 
\ Level t et al., 2006], The crossing of the equator occurs at 
13:45 local time. Two-dimensional CCD detectors span the 
wavelength range from 270 nm to 500 nm (FWHM 0.45 nm 
- 1.0 nm) in one dimension, and the other dimension is 
binned to monitor 60 different ground footprints perpendicu- 
lar to the flight direction in a push broom manner. The width 
of the resulting swath is ca. 2600 km, and global cover- 
age is obtained within one day. The size of the ground 
footprint varies across the swath from 13 x24 km 2 at 
nadir to~40x 160 km 2 for the edge of the orbit due to 
the optical aberrations and asymmetric alignment. Due 
to the so-called row anomaly caused by partially blocked 
entrance optics, some of the 60 pixels have to be 
excluded in the further analysis. 

[is] The N0 2 tropVCDs data product used here (v2.1) is 
the standard NASA OMI N0 2 product and is based on the 
algorithm described in Bucsela et al. [2013]. Briefly, the 
recorded spectra are analyzed with the DOAS method in a 
fitting window from 405 nm to 465 nm applying the 
Vandaele et al. [2002] N0 2 absorption cross section and a 
reference solar irradiance spectrum. The obtained SCDs are 
then corrected for instrumental artifacts. This is called 
destriping since the effect varies across the orbital track. To 
separate stratospheric and tropospheric components, applica- 
tion of stratospheric AMFs to destriped SCDs yield initial 
VCDs. Areas of tropospheric contamination in the 
stratospheric N0 2 field are identified using monthly mean 
tropospheric N0 2 columns from GMI simulations. Those 
regions are then masked, and the residual field of the 
stratospheric VCDs measured outside the masked regions is 
interpolated to estimate stratospheric N0 2 columns for each 
measurement. The stratospheric SCDs are subsequently 
subtracted from the original SCDs yielding the tropSCDs. 
The tropospheric AMFs (tropAMF) are calculated using a 
precomputed scattering-weight table from TOMRAD 
[Dave, 1965] and monthly mean N0 2 profiles from the 
GMI simulation. The algorithm uses OMI-based reflectivity 
[Kleipool et al., 2008], a cloud fraction and cloud pressure 
derived as described by Acarreta et al. [2004], temperature 
profiles from the GEOS-5 meteorological field, and the 
ETOP05 topography. The tropopause height is obtained 
from GEOS-5 monthly tropopause pressures. The effects of 
aerosols on the OMI N0 2 retrieval is implicitly accounted 
for through the use of the OMI-derived surface reflectivity 
which is usually larger than the true surface reflectivity due 
to scattering from aerosols [Torres et al., 2007], and through 
the OMI cloud parameters [Boersma et al., 2011]. 

[ 19 ] The uncertainty in the SCDs is ca. 10 15 molecule 
cm" 2 which corresponds to about 10% of the total slant 
column for polluted regions, and the stratospheric corrections 
leads to an uncertainty of ca. 2xl0 14 molecule cm" 2 
[Bucsela et al., 2013], Stratospheric AMFs have an 
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uncertainty of 1 - 2% and the tropAMF uncertainties range 
from ca. 20% for low cloud fraction to 30 - 80% for high 
cloud fraction, but are highly dependent on the NO 2 profile 
shape [ Wenig et al., 2008; Bucsela et al., 2013]. The overall 
error on the tropVCD is <30% under clear-sky conditions 
and typical polluted conditions (>1 x 10 15 molec cm -2 ) 
[Boersma et al., 2009; Hains et al., 2010; Irie et al., 2012; 
Bucsela et al., 2013]. 

2.3. Albedo Sensor 

[ 20 ] The University of Colorado albedo instrument con- 
sists of two four channel radiation sensors (Skye instruments 
SKR 1 850) mounted to the top and the bottom of the NOAA 
Twin Otter research aircraft pointing straight upwards and 
downwards. Each of the four telescopes is equipped with a 
custom interference filter with~10nm wide transmission 
centered at 361 nm, 479 nm, 629 nm, and 868 nm. The 
sensors are fitted with a cosine correcting diffuser plate to 
measure irradiance from the hemispherical distribution. To 
ensure the upward and downward facing channel pairs at 
any given wavelength are directly comparable, a calibration 
factor for each channel was determined by simultaneous 
zenith measurement of solar radiation during a bright sunny 
day. 1 Hz data were recorded during all flights. The ratio of 
the normalized up-welling counts to the normalized down- 
welling counts is defined as the surface albedo. 
Atmospheric backscatter for higher flight altitudes was 
corrected as follows: Radiative transfer simulations show 
that outside the aerosol layer, up-welling and down-welling 
radiation is a linear function of sensor altitude. Surface al- 
bedo measurements from high altitude flights were corrected 
for the base altitude of 1 1 00 m AGL which was assumed to 
be just outside the aerosol layer. After this correction, surface 
albedo measured during high altitude flights showed good 
agreement with low altitude flights (< 1000 m, See 
Figure 7), which were conducted to minimize the need for 
atmospheric correction. During postprocessing, the 1 Hz 
radiation data from the two sensors were averaged for 30 s 
to maximize signal to noise. The overall uncertainty on the 
surface albedo is ±5%. The instantaneous footprint of 
hemispherical irradiance measurements is a circle with 
~2.5km radius, while flying at 2.5 km, which is smaller 
than the OMI pixel size (see section 2.2). 

2.4. Flight Planning and Measurements 

[ 21 ] The Twin Otter flight plan included usually a morning 
flight (ca. 10 - 14 h local Pacific Daytime Time) with 
emphasis on a combination of low approaches to obtain pro- 
file information and constant flight altitude at ca. 2 - 2.5 km 
AGL, just outside the boundary layer, followed by an 
afternoon flight (ca. 15 - 19 h) with a constant flight altitude 
of~3.5 - 4 km AGL. Overall, we performed 51 research 
flights on 30 days over the two month period, but only 41 
flights took place on the same days as OMI measurements. 
This is due to the row anomaly of this instrument (see 
section 2.2). Coordination with satellite overpasses was not 
the primary focus of the flight planning. Rather, the frequent 
OMI overpasses coupled with the large swath width of the 
OMI instrument result in a high number of flights that are 
suitable for satellite validation. 


3. Satellite Validation 

3.1. Comparison of All Data 

[ 22 ] In a first step, tropVCDs are grouped to calculate the 
median values for specific regions over California covering 
the whole two months, i.e., individual measurements of the 
two observational platfonns are not directly compared, but 
rather the respective medians of distributions within a 
specific region. The regions, as outlined by boxes in 
Figure 2, are chosen at locations with sufficient data points 
of AMAX-DOAS measurements. These mainly urban areas 
include the Greater Los Angeles Area (boxes 1 - 5), 
Bakersfield (boxes 6 - 7), and Sacramento (boxes 8 - 10). 
OMI tropVCDs are selected when the center of a pixel is 
located within a box on the same day a Twin Otter flight took 
place over this box. Days on which either data set is not 
available are excluded. There is no other temporal coinci- 
dence criterion applied resulting in possible time differences 
between the AMAX-DOAS measurements of up to 7h 
before and 5h after the satellite overpass. Figure 3 shows 
the medians split by region for the two instruments. Both data 
sets follow a similar trend: high NO 2 tropVCDs are observed 
over the Los Angeles Metropolitan Area (boxes 1 and 2) 
where many major motorways intersect. The Los Angeles 
International Airport (LAX) and the Port of Los Angeles 
are both situated in box 2. The Inland Empire also displays 
very high N0 2 (box 3). Lowest N0 2 is found west of 
Bakersfield (box 7). In general, the regions around 
Bakersfield and Sacramento are cleaner than the Greater 
Los Angeles Area. As expected, NCL amounts are closely 
related to population and transportation. The population of 
the Greater Los Angeles Area is 17,877,006 (2010 US 
census, http://2010.census.gov/2010census/data/). In com- 
parison, the population of Sacramento is 466,488 and 
Bakersfield 347,483. In the so-called South Coast Air Basin 
(SCAB), a geopolitical area defined for the purpose of air 
quality management and delimited by the Pacific Ocean to 
the west and mountain ranges in all other directions, on-road 
motor vehicles and other mobile sources together account for 
91% of the total NO x emissions (Emission data for 2008 by 
California Air Resources Board: http://www.arb.ca.gov/ 
app/emsinv/fcemssumcat2009 .php) . 

[ 23 ] The satellite data are lower than the AMAX-DOAS 
measurements for enhanced NO 2 , see boxes 1 and 3. This 
is also reflected in the slope of 1 .4 ± 0. 1 for the linear regres- 
sion analysis of the two data sets (Figure 4, CU-AMAX- 
DOAS vs. NASA OMI, Pearson correlation coefficient 
R=0.97, offset — 0.7±0.5 x 10 15 molecule cm -2 ). If the data 
from boxes 1 and 3 are omitted, the slope is reduced to 
0.9 ±0.1 and hence insignificantly different from 1. The 
Pearson correlation coefficient is 0.96 and the offset 
0.4 ± 0.3 x 10 15 molecule cm -2 . Please note, the simple linear 
regression does not take the uncertainty of the individual data 
points into account, and these offsets are much smaller than 
the range of the 25 th and 75 th percentile of the data (see 
Figure 3) or the uncertainty of the stratospheric correction. 
The reason for the underestimation of the OMI data for 
enhanced N0 2 most likely originates from the selection 
criteria for the satellite pixels. As stated above, the center of 
the ground footprint has to be within the box. Furthermore, 
also pixels at the side of the OMI swath are included and 
those can be up to 160 km wide. This results in a bias to 
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Figure 2. Overview of CU-AMAX-DOAS NO 2 tropVCDs over California 19 May - 19 July 2010. To 
show the variation in the background NO 2 , the color scale was set to a range up to 1 x 10 16 molecule 
cm“ 2 , but highest observed NO 2 was ~3 x 10 16 molecule cm~ 2 . White dashed boxes in panel a mark the 
outlines of panels b-d. The dotted boxes (1-10) in b-d outline the areas for which the data were averaged 
for comparison with NASA OMI NO 2 averages as summarized in Figures 3 and 4. Stars in panel b mark the 


location of profiles as presented in Figure 8; from 
Banning airports. 

background values since the actual area sampled by OMI 
included in the average will extend beyond the delimiters 
of the boxes by up to 80 km. Especially, the highly polluted 
boxes are surrounded by cleaner areas, e.g., box 3 is adjacent 
to the San Gabriel Mountains and the High Desert located to 
the north and west. 

[ 24 ] Only for the Greater Los Angeles Area (boxes 1 - 5), 
statistics were good enough to separate the data into weekend 
and weekday measurements. A clear decrease can be observed 
(see Figure 3) which is caused by a reduction of heavy-duty 
diesel-fueled vehicles and is consistent with previous studies 
for California [e.g., Marr and Harley, 2002; Chinkin el al., 


west to east: Santa Monica, Brackett, Ontario, and 


2003; Pollack et al., 2012 and references therein]: the NO 2 
is lower by 38±24% and 33± 11% for AMAX-DOAS and 
OMI, respectively. In the regression analysis (see Figure 4), 
the slope for the weekday data of 1.8 ±0.3 is clearly larger 
than the 1: 1-line, whereas the weekend regression with a 
slope of 1.1 ±0.2 is consistent with a slope of unity. Both 
offsets are insignificantly different from zero. However, the 
two slopes are statistically different indicating as already 
mentioned above, the overestimation of the AMAX-DOAS 
data in comparison to the OMI measurements is nonlinear. 
The reason for this is the exponential increase of the satellite 
pixels towards the edge of the swath. 
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Figure 3. Comparison of NASA OMI and CU-AMAX-DOAS median tropVCDs. Data were averaged 
over two months and within the boxes 1-10 as indicated in Figures 2b-d. Boxes 1-5 are located over 
the Greater Los Angeles Area, 6-7 over Bakersfield, and 8-10 over Sacramento. A linear regression 
analysis was performed on the median of the data, and the results are shown in Figure 4. 
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N0 2 tropVCD [10 1S molec cm' 1 ] 

Figure 4. Linear regression analysis for data from F igure 3 . 
The correlation coefficient R and the slopes and offsets are 
given in the legend. The offset has to be multiplied by 
[10 15 molecule cm -2 ]. 

3.2. Comparison Based on Individual OMI Pixels 

[25] The large number of data points of the CU-AMAX- 
DOAS instrument allows for a validation of the NASA 
OMI data on a pixel level. Figure 5 shows an example flight 
on 1 June 2010 where the AMAX-DOAS tropVCDs are 
plotted on top of the gridded OMI data. The individual 
AMAX-DOAS measurements can be combined within the 
area of a satellite pixel to calculate a mean AMAX-DOAS 
NO2 tropVCD which can then be compared to the specific 


OMI tropVCD. In order to identify the AMAX-DOAS data 
points within that area that is given by the four comer points 
of the OMI pixel, the so-called point-in-polygon problem 
was approached with a ray-casting algorithm. The flight 
presented in Figure 5 is from earlier in the campaign where 
the detector exposure times were 20 s. Nevertheless, there 
are up to four AMAX-DOAS measurements per OMI pixel. 

[26] Several factors can affect the quality of the compari- 
son exercise. In the following, we discuss five key parame- 
ters: (1) Since the NO2 lifetime is in the range of a few 
hours and also transport can take place, the time difference 
between the measurements is very important. (2) Although 
most sources of NO2 are close to the surface, our vertical 
profile measurements during this campaign have shown that 
the NO2 can be mixed in the BL up to ca. 2 km [Baidar 
et ah, 2013] and hence the aircraft altitude determines the 
fraction of the NO2 tropospheric column sampled by the 
AMAX-DOAS. (3) As was observed above, the quality of 
the comparison suffers when the large pixels at the sides of 
the swath are included, and hence this parameter is tested 
here as well. (4) Obviously, the cloud radiance fraction 
determines how representative the measured OMI column 
is for the whole pixel. In general, the cloud radiance fraction 
during the campaign was rather low though. 

[27] (5) Further, we define a quantity we call the normal- 
ized distance to assess the representativeness of the 
AMAX-DOAS measurements within a pixel with respect to 
the sampled area. The OMI pixels are divided into similarly 
sized subgrid cells with a side length of~7.5km each 
resulting in 3x2 to 20 x 4 dimensional subgrids. This 
variation is due to the increasing footprint of the OMI 
measurements from the middle to the side of the swath. In a 
first step, it is checked whether the above defined subgrid 
cells contain any data points. Matrices with the dimension 
of the subgrid are defined for each of these individual cells. 



0.1 1 2 3 A 5 6 7 9 10 


Figure 5. CU-AMAX-DOAS N0 2 vertical columns plotted on top of a N0 2 NASA OMI swath for a 
selected flight on 1 June 2010. The white plus indicates the AMAX-DOAS measurement recorded at the 
time of the OMI measurement, circles indicate a time difference <30 min, diamonds >30 min, and the size 
of the symbols is inversely proportional to the time difference. 
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Figure 6. Panel f summarizes the data from Figure 5. The letter and number combination of the x-axis 
indicates the location of the individual OMI pixel in the grid of Figure 5. Panel a: normalized distance 
(see text and Figure Al), panel b: cloud radiance fraction of the OMI measurement, panel c: mean time 
difference between the AMAX-DOAS and the OMI measurement, panel d: Twin Otter altitude, panel e: 
number of individual AMAX-DOAS data points that were in the average of panel f. The red lines and 
associated arrows (green: good, red: less optimal) indicate chosen thresholds for the selection criteria, 
and the blue dots below panel f mark the good data points based on the selection criterion which contributed 
to the regression analysis shown in Table 2. 


The step distance between occupied subgrid cells are 
assigned as the matrix elements for the individual occupied 
cell matrix. For example, the distance between adjacent 
occupied cells is 1, between diagonally adjacent occupied 
cells 2, a knight's move distance is 3, etc. These individual el- 
ements are summed up and divided by the value of the best 
case scenario, i.e., every subgrid cell contains at least one 
data point. Consequently, the normalized distance is a 
number between 0 and 1 and an approximate measure for 
the distribution of the AMAX-DOAS data points within an 
OMI pixel. How to obtain the normalized distance is 
sketched in Figure Al for an example of three occupied cells 
of a 3 x 2 subgrid. This normalized distance is less dependent 
on the actual number of data points which varied during the 
campaign due to different integration times. 

[ 28 ] The data for the example flight of Figure 5 are summa- 
rized in Figure 6f with each spatially coincident measure- 
ment identified by a letter and number coordinate as 
defined in Figure 5. Figures 6a - e show the key parameters 
as described above: The normalized distance of the 
AMAX-DOAS data points (Figure 6a), the OMI cloud 


radiation fraction (Figure 6b), the mean time difference 
between the AMAX-DOAS measurements and the OMI 
measurement (Figure 6c), the aircraft altitude (Figure 6d), 
and the number of AMAX-DOAS data points (Figure 6e). 
During this example flight, the included data points are from 
the OMI pixel numbers 18 to 24 (with 30 and 31 being the 
center of the orbit) within the OMI swath. With a combina- 
tion of selection criteria as indicated by the horizontal red 
lines in Figures 6a - d (nonnalized distance >0.01, cloud 
radiance fraction <50%, OMI overpass within ±2 h, flight 
altitude >2 km), the initial 31 data points are reduced to 13 
as indicated by the blue marks below panel f. The resulting 
slope of the linear correlation is 0.8 ±0.2, with a moderately 
strong correlation coefficient of R = 0.78 and an offset of 
1.1 ±0.9 x 10 15 molecule cm~ 2 . Please note the uncertainties 
of the measurements are neglected in the regression since the 
uncertainty for each individual data point is unknown (see 
sections 2.1 and 2.2). This and all following regression 
results are based on the AMAX-DOAS data acting as 
ordinate and the OMI data as abscissa in the Cartesian 
coordinate system. 
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Table 1. Regression Analysis for All Data on the Pixel Level a 


Scenario 1 * 

Slope 

Offset 

[10 15 molec/cm 2 ] 

N OMI Pixels 

(AMAX-DOAS Data Points) 

R 

Standard Deviation of Fit d 
[10 15 molec/cm 2 ] 

All data 

1.12 ±0.04 

0.67±0.13 

1016(34,058) 

0.69 

2.9 

Pixels 

6-55 

1.07 ±0.04 

0.25±0.14 

546 (10,555) 

0.78 

2.4 

11-50 

1.04 ±0.04 

0.35±0.19 

363 (5300) 

0.78 

2.6 

16^15 

1.06 ±0.05 

0.37±0.20 

284 (3507) 

0.79 

2.5 

Normalized distance : 

>0.01 1.17±0.04 

0.36±0.15 

602 (28,559) 

0.78 

2.5 

>0.05 

1.21 ±0.05 

0.23 ±0.24 

312 (20,649) 

0.80 

2.8 

>0.1 

1.30 ±0.08 

— 0.09±0.35 

187 (14,756) 

0.79 

3.0 

Cloud radiance fraction . 

<50% 1.14±0.04 

0.51 ±0.14 

908 (32,153) 

0.70 

2.8 

<30% 

1.14 ±0.04 

0.51 ±0.11 

848 (28,135) 

0.70 

2.8 

OMI overpass'. 

±2h 

1 . 1 8 ± 0.07 

0.69±0.23 

382 (12,607) 

0.66 

3.1 

±1 h 

1 .30 ± 0. 12 

0.23±0.36 

131 (3248) 

0.69 

2.8 

Flight altitude : 

>2 km 

1.05 ±0.04 

0.58 ± 0. 13 

830 (27,994) 

0.69 

2.5 

>3 km 

1.12 ±0.06 

0.31 ±0.20 

391 (12,208) 

0.72 

2.8 

>4 km 

1.39 ±0.08 

— 0.27±0.26 

126 (4177) 

0.84 

2.0 


“AMAX-DOAS data acting as ordinate and the OMI data as abscissa in a Cartesian coordinate system. 

“See section 3 for details. 

“Numbers refer to the pixel number i n the OMI swath. 

d The standard deviation of the fit = ^/ v [ , f'\y t — + ^4)] 2 with N being the number of measurements, x, and y, the measurements themselves, and A and 

B the intercept and slope of the regression, respectively. 


[ 29 ] In the following, this kind of regression analysis is 
applied to all of the available data. As before, the initial selec- 
tion criterion is that the NASA OMI and the CU- AMAX- 
DOAS measurements were performed on the same day. 


This results in 1016 individual OMI measurements being 
initially included in the comparison and 34,058 colocated 
AMAX-DOAS measurements for the two months of the 
measurement period. For the comparison of the two data sets, 


Table 2. Regression Analysis Results for All Data on the Pixel Level for a Combination of the Best-Off - Scenarios Based on Table l a 


Scenario b 

Slope 

Offset 

[10 15 molec/cm 2 ] 

N OMI Pixel 
(AMAX-DOAS Points) 

R 

Standard Deviation of Fit d 
[10 15 molec/cm 2 ] 

All data 

1 . 12 ± 0.04 

0.67±0.13 

1016(34,058) 

0.69 

2.9 

Pixels 1 1-50“ 

1.04±0.04 

0.35±0.19 

363 (5300) 

0.78 

2.6 

Pixels 1 1— 50 c 
+ Nonn. distance >0.01 

1.03 ±0.04 

0.35±0.16 

282 (5021) 

0.86 

2.0 

Pixels 1 1-50“ 

+ Nonn. distance >0.01 
+ Cloud rad. fraction <50% 

1.01 ±0.04 

0.30±0.17 

252 (4671) 

0.86 

2.0 

Pixels 1 1-50“ 

+ Nonn. distance >0.01 
+ Cloud rad. fraction <50% 
+ OMI overpass ± 2 h 

1.06±0.08 

0.20±0.37 

81 (1220) 

0.82 

2.3 

Pixels 1 1— 50 c 
+ Nonn. distance >0.01 
+ Cloud rad. fraction <50% 
+ OMI overpass ± 2 h 
+ Flight altitude >2 km 

1.03±0.09 

0.14±0.43 

56 (632) 

0.85 

2.3 


“AMAX-DOAS data acting as ordinate and the OMI data as abscissa in a Cartesian coordinate system. 
b See section 3 for details. 

“Numbers refer to the pixel number i n the OMI swath. 

d The standard deviation of the fit = ^ ~y t — ( Sx, + A)] 2 with N being the number of measurements, Xj and y, the measurements themselves, and A and 
B the intercept and slope of the regression, respectively. 
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Figure 7. Surface albedo at 479 nm as measured by the radiometers on board the NOAA Twin Otter on 3 1 
May 2010 (symbol: diamonds) and 12 July 2010 (symbol: discs) in the Los Angeles region. 


the AMAX-DOAS measurements are first averaged within 
the subgrid cells as defined for the calculation of the normal- 
ized distance, and the resulting mean values are then again 
averaged. This two-step calculation of the mean is done to 
avoid any weighting towards heavily sampled areas of the 
OMI pixel. Table 1 shows the results for the linear regression 
analysis for the averaged AMAX-DOAS measurements 
compared to the OMI data and how the results change when 
constraining the above described key parameters. 

[ 30 ] In general, the correlation coefficient increases with 
more stringent limits and so does the slope up to values 
significantly higher than 1 (see Table 1), i.e., the AMAX- 
DOAS measurements are higher than the NASA OMI ones 
(by up to 39 ±8%). The only exception to the increase of 
the slope with stricter coincidence criteria is when removing 
the large OMI pixels at the sides of the swaths from the 
comparison; then the slope actually decreases (although not 
significantly within the uncertainty range). The offset does 
not seem to follow a certain pattern, but is mostly positive. 
However, the numbers are smaller than the uncertainty on 
the stratospheric correction of 0.5 - 2 x 10 15 molecule cm -2 
for the CU-AMAX-DOAS measurements and 0.2 x 10 15 
molecule cm -2 for the NASA OMI. This somewhat 
surprising trend in the slope can be explained when viewed 
in relation to the results described in the following exercise. 

[ 31 ] Table 2 and Figure A2 summarize the effect of succes- 
sively applying the weakest selection criteria to the full data 
set of the initial 1016 OMI pixels again. The weakest criteria 
were chosen to leave sufficient data points for good statistics. 
After removing the sides of the swath, adding any other 
limitation results in a strong correlation of R> 0.82 and does 
not change the slope significantly from 1 . In the end, for the 
remaining pixels, the agreement is excellent, and the 
concluding slope is 1.03 ±0.09. After removing the large 
pixels of the side of the swath, constraining the nonnalized 
distance decreases the standard deviation of the fit and 
increases the correlation coefficient showing that the normal- 
ized distance is indeed an important factor. Limiting the cloud 
radiance fraction only results in negligible effects since the 
pixels had more or less already a very low cloud fraction to 
start with. Reducing the difference in the OMI overpass time 
to ±2h within the AMAX-DOAS measurements removed 


many of the data points, and the results are not very conclusive 
anymore. The standard deviation of the fit as well as the uncer- 
tainty on the slope significantly increased. The consecutive 
removal of measurements with an aircraft altitude <2 km does 
not change the picture anymore which suggests that most of 
the NOt is below 1.8 km, the initial selection criterion. 

[ 32 ] As mentioned above, these results seem to be in 
contrast to Table 1 where the individual limiting criteria 
mostly increased the slope, i.e., the AMAX-DOAS data seem 
to overestimate the NO 2 in comparison to the OMI pixel. 
This previous regression analysis was performed including 
all pixels of the OMI swath. The area of the footprint of a 
nadir AMAX-DOAS pixel is 0.011km 2 when flying at 
4 km altitude and integrating for 2 s (see section 2.1 last 
paragraph). This is relatively small in comparison to the 
OMI footprint area of 3 12-6400 km 2 . This indicates that in 
combination with the small FOV of the AMAX-DOAS 
instrument, the NOAA Twin Otter is not fast enough to 



0 4 8 12 0.00 0.05 0.10 0.15 0.20 


N0 2 mixing ratio [ppb] aerosol extinction coeff. [km 1 ] 

Figure 8. Profiles for NO 2 (panel a) and the aerosol 
extinction coefficient at 477 nm (panel b) as derived from 
CU-AMAX-DOAS measurements on 16 July 2010 [Baidar 
et al, 2013]. The locations of these profiles are marked in 
Figure 2b, and they span the SCAB in roughly east-west 
direction from the Banning mountain pass to the coast 
(Santa Monica Airport). The profiles were normalized to 
0 m ASL. The red line shows the N0 2 GMI model profile 
for Los Angeles for July 2010 multiplied by factor 5. 
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Table 3. Satellite tropAMF Sensitivity Studies for 17° SZA and 435 nm a 


Scenario 





N0 2 Profile 




Banning 

Model 

Brackett 

Model 

Ontario 

Model 

Santa 

Monica 

Model 

0 a 

Measurements 

0 a 

Model 

Albedo 0.083, no aerosols 

1.28 

1.36 b 

1.21 

1.36 

1.12 

1.36 

1.10 

1.36 

0 a 1.18 

1.36 b 

Albedo 0.1, aerosol extinction from Figure 8b 

I.38 c 

1.45 

1.32 c 

1.46 

1.33 c 

1.46 

1.21° 

1.45 

0 a 1.31° 

0 a 1 46 

Albedo 0. 1 , no aerosols d 

1.37 

1.44 

1.32 

1.44 

1.23 

1.44 

1.22 

1.44 

0 a 1.29 

1.44 


a 0: mean. 

b tropAMF calculated with the same N0 2 profile and same settings for the aerosol/surface albedo parameterization as the actual OMI tropAMF calculations. 
TropAMF calculated using the Twin-Otter-measured average aerosol extinction, average N0 2 profiles (Figures 8a and 8b), and albedo. 
d Scenario to illustrate that aerosols have only little influence. 


representatively sample the large OMI pixels in a sufficiently 
short time interval. Since the flight tracks were mainly over 
polluted areas, the MAX-DOAS measurements are automat- 
ically biased to the enhanced NO 2 . In order to verify this, 
even more satellite data are needed so that the regression 
presented in Table 1 could be amended to only include 
OMI data with pixel numbers of 6-55 or even 11-50. 
However, OMI pixels in the same swath and adjacent swaths 
are measured almost simultaneously. Hence, to improve on 
this validation technique, extended time periods of measure- 
ments, several aircrafts measuring simultaneously, a larger 
FOV of the AMAX-DOAS instrument, faster aircrafts, or a 
combination of the four is needed. 


4. Sensitivity Studies on the Satellite AMF 

[ 33 ] hi this section, the assumptions of the satellite radiative 
transfer are tested with available auxiliary measurements of 
trace gas profiles, aerosol extinction profiles, and surface 
albedo. As can be seen from Figure 1 , the sensitivity of the mea- 
surements towards the actual NO 2 profile changes with altitude. 
Hence, the tropAMF depends on the a priori trace gas profile: 


tropAMF 


'ZiBAMFi x VCD a i priori 
tropVCD" priori 


( 2 ) 


with BAMFj and VC IF/ priori being the BAMF and the trace 
gas a priori partial vertical column of the layer i and 
tropVCD a pnon the integrated a priori profile. The summation 
is performed from the surface to the tropopause for discrete 
layers. If the true profile is different from the a priori profile, 
this can introduce significant errors in the tropVCD retrieved 
with this tropAMF. 

[ 34 ] As mentioned above, the satellite retrieval uncertainty 
is dominated by the uncertainty in the tropAMF calculation 
and there the most important parameters are trace gas con- 
centration profile and surface albedo [Boersma et al., 2004; 
Richter and Burrows, 2002], By using independent profile 
observations measured by the AMAX-DOAS instrument 
while flying low approaches over airports [see Baidar 
et al., 2013] as well as surface albedo by the radiometers, 
sensitivity studies were perfonned to test the assumptions 
made in the NASA OMI retrievals. As mentioned above, 
the radiative transfer calculation of the OMI retrievals param- 
eterizes the aerosol load in the atmosphere in combination 
with the surface albedo by applying the so-called effective al- 
bedo calculated from the measurement itself. Here, a case 
was chosen for an OMI nadir pixel for Los Angeles for 
July 2010 which was deemed as relatively low cloud by the 


OMI retrieval algorithm having an effective albedo of 
0.083. The corresponding N0 2 profile is given in Figure 8a 
and was obtained from the GSFC GMI CTM model. The 
resulting tropAMF is 1.41 at 435 nm. 

[ 35 ] Figure 7 shows the surface albedo as measured from 
the Twin Otter in summer 2010 over the Greater Los 
Angeles Area, parts of the High Desert and over the San 
Bernardino Mountains. Mostly, the surface albedo is around 
0.10 ±0.02. Only over some of the mountain ranges and the 
Pacific Ocean, the surface albedo is lower: ca. 0.03 - 0.07. 
Figures 8a and 8b show the NO 2 and the 477 nm aerosol 
extinction profiles measured during four low approaches 
over the Greater Los Angeles Area (location marked in 
Figure 2b) on 16 July 2010. Both shape and absolute values 
of the trace gas profiles differ significantly from the 
model profile. 

[36] Here, tropAMFs are calculated for the OMI nadir 
viewing geometry at a SZA of 17° for a clear-sky scenario. 
Simulations are perfonned for the five different NO 2 profiles 
from Figure 8a. The McArtim code is employed because 
aerosol extinction can be treated explicitly. The aerosol 
optical properties were chosen as 0.94 for the single scatter- 
ing albedo, typical for wavelengths >400 nm, and 0.68 as 
asymmetry parameter, typical for polluted environments. 
These are the setting used for retrieving the NCB and aerosol 
profiles shown in Figure 8 using the same radiative transfer 
model McArtim. The aerosol extinction profiles correspond- 
ing to the NO 2 profiles at the four different locations are 
applied (see Figure 8b). Calculations were also performed 
for a Rayleigh atmosphere with 0.083 (corresponds to the 
effective albedo) and 0.1 surface albedo. Profiles for 
pressure, temperature, and ozone were taken from the US 
standard atmosphere to be consistent with what was used 
during the profile retrievals [Baidar et al., 2013]. However, 
using the GEOS-5 meteorological field (which drives the 
GMI model) pressure and temperature profiles increases the 
tropAMFs by ~1 — 2%. The tropAMF obtained by mimicking 
the OMI retrievals by using the albedo settings (i.e., effective 
albedo of 0.083 and no aerosols) and the NO 2 model profile 
results in a tropAMF of 1 .36 which is within 4% of the actual 
OMI tropAMF applied to the measured slant columns. The 
aim of these radiative transfer calculation comparisons is 
not to reproduce the exact OMI tropAMF, but rather to study 
the differences to the standard scenario, i.e., 1.36 tropAMF in 
order to identify the main contributors as well as obtaining an 
estimate for the satellite tropAMF uncertainty. 

[ 37 ] Table 3 summarizes the results for the different 
scenarios. On average, the GMI model profile results in a 
larger tropAMF in comparison to the measured NO 2 profile 
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tropAMF by 11-15%. An increase in the surface albedo of 
20% results in an average tropAMF increase of 9% for the 
measured profiles. When using the GMI model profile in 
combination with the different measured aerosol extinction 
profiles, the tropAMF changes by 51%. The weak depen- 
dence of the satellite (but also the AMAX-DOAS) 
tropAMF on the aerosol extinction profile was already 
explained in the introduction. Using a combination of the 
Twin Otter measurements for surface albedo, aerosol 
extinction profile, and N0 2 profile and giving equal weights 
to the four individual profile combinations yields on average 
a tropAMF of 1.31. This value is about 4% lower than the 
tropAMF of 1.36 using the OMI settings, i.e., an effective 
albedo of 0.083 in combination with a Rayleigh atmosphere, 
and the model N0 2 profile. Looking at the individual 
measured profiles, the maximum deviation is for the coastal 
profile at Santa Monica with 12%. It is worth noting that this 
value of 12% is smaller than the differences caused by only 
using the different N0 2 profiles, i.e., up to 15%. Applying 
the measured effective albedo seems to partly offset the 
nonideal choice of the model trace gas profile. 

[38] In summary, the NASA OMI tropAMF for the Los 
Angeles area for July 2010 is only as good as the estimate 
for the albedo and especially the trace gas profile. 
Considering the assumptions in the simulations and the range 
of values, a tripling of the average difference seems to be 
appropriate to estimate the satellite tropAMF uncertainty 
for the Los Angeles area for the summer to be 12%. This is 
much better than the estimate described in section 2.2 for 
OMI retrievals in general (ca. 20% for low cloud fraction, 
30 - 80% for high cloud fraction [Wenig et al., 2008; 
Bucsela et al., 2013]), but is consistent with the results of 
the linear regression study in the previous section. 

5. Additional Error Sources 

[ 39 ] Error sources in the OMI and the AMAX-DOAS data 
products which have not been discussed so far in this 
manuscript are the effect of the orography on the field in 
the RTM calculations, the temperature dependence in the 
N0 2 absorption cross section, and the contribution to the 
N0 2 column from above the aircraft altitude: The radiative 
transfer models used in this study treat the atmosphere and 
the surface elevation as isotropic. However, the terrain height 
of the SCAB is highly variable. On the other hand, areas with 
highest N0 2 are usually not in the foothills of the surround- 
ing mountains. Also, this effect should become significant 
only for large SZAs in combination with the position of the 
sun in the direction of a mountain range. Hence, the influence 
is most likely negligible, especially when investigating statis- 
tical ensembles of data. A detailed study is beyond the scope 
of this manuscript. 

[ 40 ] The AMAX-DOAS analysis was performed with a 
combination of two N0 2 cross sections at 220 K and 294 K. 
When using only the cross section at 294 K in the analysis, 
the differences in those dSCDs to the dSCDs obtained with 
the combination of the two absorption cross sections are neg- 
ligible. However, using the 220 K cross section yields 20% 
smaller dSCDs. This highlights the need for treating the 
temperature dependence of the N0 2 absorption cross section 
in the two retrievals in a similar way for a meaningful 
comparison. In the satellite algorithm, the temperature 


dependence of the cross section is treated explicitly by assim- 
ilating a temperature profile and calculating a correction 
factor for the cross section [Bucsela et al., 2013]. The 
average temperature in the lowest kilometer in the model 
profile for July over California is 295 K and therefore close 
to the value in the AMAX-DOAS retrievals. 

[ 41 ] N0 2 can be produced in the free troposphere from 
lightning, but California was mainly cloud free during the 
campaign. However, the GMI model profile yields a partial 
column of 4 x 10 14 molecule cm -2 from 5 km to the tropo- 
pause. This partial column above the aircraft is comparable 
to the offsets found in Table 1. 

[ 42 ] The CU- AM AX-DO AS retrieval error is already low, 
but could possibly be improved for higher SZA by using op- 
timal estimation to obtain a vertical column including obser- 
vations from additional viewing angles. This will be 
examined in a follow-up research paper. 


6. Summary and Conclusions 

[ 43 ] In this paper we presented (1) a comparison of NASA 
OMI with CU-AMAX-DOAS N0 2 median tropVCDs over 
2 months for defined areas and a comparison on an individual 
OMI pixel basis over California in summer 2010, and (2) 
sensitivity studies on the satellite AMF. The results can be 
summarized as follows: 

[ 44 ] 1. Fair agreement was found for the temporally and 
spatially averaged data comparison (section 3.1). The slope 
between the tropVCDs of the two instruments is 1.4 ±0.1 
and the correlation coefficient 0.97. This is caused by a 
combination of large OMI pixels at the side of the swath 
and the SCAB being surrounded by relatively unpolluted 
areas. Therefore, large pixels should be excluded if the 
coincidence criterion is based on pixel center coordinates. 

[ 45 ] 2. The mean AMAX-DOAS tropVCDs were calcu- 
lated for coincident OMI pixels. A regression analysis was 
performed and the successive application of a combination 
of individual selection criteria to the data (i.e., pixel number 
in the OMI swath, normalized distance, cloud radiance frac- 
tion, satellite overpass time, and aircraft altitude) led to a 
strong correlation of 0.85 and a slope of 1.03 ±0.09 in the 
end showing good agreement. The main driver seems to be 
the removal of the large pixels from the comparison. The 
AMAX-DOAS footprint area of -0.01 lkm 2 is relatively 
small in comparison to the OMI footprint area of 3 12 — 
6400 km 2 . Optimizing the normalized distance is also impor- 
tant. The cloud radiance fraction and the aircraft altitude have 
only minor impacts on the results since the former is very low 
anyways over Los Angeles in summer and in the latter case 
since most of the tropospheric N0 2 seems to be confined 
below 1800 m. 

[46] 3. The statistics of this validation technique can be 
improved by extending the time periods of measurements, 
coordinating several aircrafts at the same time, developing 
an AMAX-DOAS instrument with a larger FOV, deploying 
a faster aircraft or a combination of the four. 

[ 47 ] 4. Sensitivity studies on the satellite AMF showed that 
the radiative transfer is rather independent of the range of 
aerosol load as encountered during the campaign, but highly 
dependent on the surface albedo and the trace gas a priori 
profile shape. 
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[48] 5. The uncertainty of the NASA OMI tropAMF is 
estimated to be 12% for summer over the Los Angeles area. 
Notably, the area probed is characterized by a rather high 
surface albedo (here 10% at 479 nm) and low AOD. A 
generalization of the satellite uncertainty over areas with a 
different surface albedo may not be straightforward. 

[ 49 ] 6. The observed weekly cycle of the NO 2 is consistent 
with previous studies. A decrease of 38 ± 24% and 33 ± 1 1% 
for AMAX-DOAS and OMI measurements, respectively, is 
found for the weekend in comparison to weekdays. While 
the effects on ozone are not the subject of this study, we note 
that this decrease in NO 2 leads to higher ozone during week- 
ends in June (consistent with Pollack et al. [2012], and ref- 
erences therein), but lower ozone during a hot weekend 
case in July [Baidar et al., 2012], 

[ 50 ] The aircraft observations reported here were only par- 
tially intended for satellite validation. Purposeful validation, 
flying grids coincident in time and space with the satellite 
overpass, and adjustments to the actual AMAX-DOAS 
instrument could provide a more detailed data set. All this 
said, it can be concluded that (1) CU-AMAX-DOAS is well 
suited for satellite validation when the above precautions 
are taken into account, and the data collected over 
California are a valuable data set to validate other satellite 
instruments and (2) the NASA OMI tropospheric N0 2 
product (v2.1) delivers high quality data for the Californian 
summer season. 
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